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ABSTRACT 

We present the results of a search for optical counterparts to the two quiescent low mass X-ray binaries (X5 
and X7) in the globular cluster 47 Tucanae, using high quality Chandra and //^T images. A faint blue (V = 21.7; 
U -V = 0.9) star within Of'03 of the eclipsing system X5 shows variability on both short and long timescales, and 
is the counterpart of the X-ray source. The colors and variability of this object are consistent with the combination 
of light from an accretion disk and a red main sequence star (possibly somewhat larger than a normal MS star with 
similar luminosity). No evidence is found for a star showing either variability or unusual colors near the position 
of X7, but a probable chance superposition of a star with V = 20.25 limits the depth of our search. 

Subject headings: binaries: general - globular clusters: individual (47 Tucanae) - techniques: photometric - 
X-rays: binaries 



I 1. INTRODUCTION 

It has long been thought that quiescent low mass X-ray bina- 
ries (qLMXBs) dominate the most luminous of the dim sources 
in globular clusters (Hertz & Grindlay 1983, Verbunt et al. 
^ 1984). Recent observations using Chandra/ACIS imaging and 
. spectroscopy have demonstrated that this is indeed the case. 
Two systems, X5 and X7 in the massive globular cluster 47 Tuc 
were previously suspected to be qLMXBs (Hasinger, Johnston 
& Verbunt 1994, Verbunt & Hasinger 1998) but the sensitivity 
and resolution of Chandra was required to confirm this suspi- 
cion (Grindlay et al. 2001a, hereafter GHEOla and Heinke et 
al. 2001a, in preparation, hereafter HGLOl). One qLMXB has 
also been found in each of NGC 6397 (Grindlay et al. 2001b; 
hereafter GHEOlb) and co Cen (Rutledge et al. 2001). These 4 
qLMXB systems all have thermal spectra that are well modeled 
by hydrogen atmospheres of hot neutron stars (NSs), with no 
power law components required. None of them are obviously 
variable with the exception of X5 (HGLOl) which shows deep 
eclipses as well as dips showing increased neutral hydrogen (X7 
shows marginal evidence for a 5.5 hr period). 

The logical extension of this work is to search for optical 
counterparts to these sources, using the potent combination of 
Chandra and HST. With astrometric errors < Of'l routinely be- 
ing achieved for X-ray sources, optical identifications are being 
reported with unprecedented frequency. These identifications 
include cataclysmic variables (CVs) and BY Draconis variables 
(GHEOla, GHEOlb), milHsecond pulsars (MSPs; Edmonds et 
al. 2001, hereafter EGHOl and Ferraro et al. 2001) and active 
LMXBs (Heinke et al. 2001b and White & Angelini 2001). 

Searches for optical counterparts to qLMXBs have been less 
successful. No counterpart has been found for the NGC 6397 
qLMXB (any possible companion has My> 1 1 ; GHEOlb) while 
the CO Cen qLMXB lies outside the field of view (FOV) of cur- 
rent HST datasets and stellar crowding will limit deep searches 
from the ground. Here, we report the use of high quality Chan- 
dra and HST data to search for optical counterparts to the 47 
Tuc qLMXBs X5 and X7. We have discovered a faint, blue and 
variable counterpart to the eclipsing X5, as reported briefly in 
GHEOla. This detection, combined with the well determined 
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period, distance, inclination and X-ray spectrum of X5 makes 
this the best constrained qLMXB known. We also report limits 
on the qLMXB X7. The astrometry, photometry and time series 
for both of these searches are described below. 

2. OBSERVATIONS AND ANALYSIS 

Details of the Chandra data used here are given by GHEOla 
and HGLOl. The qLMXBs X5 and X7 have 4576 and 5488 
counts respectively (over the 72 ksec observation), with inter- 
nal, Icr errors of 0f'0082 and Of'0089 respectively. To search 
for optical counterparts to X5 and X7, two HST datasets have 
been analysed, the 8.3 d observations of Gilliland et al. 2000 
(GO-8267: July 3 1999 to July 11 1999) and the archival data 
of Meylan obtained in three different epochs with ~ 2 year 
spacings (GO-5912: October 25 1995; GO-6467: November 
3 1997; GO-7503: October 28 1999). The Gilliland data pro- 
vides exquisite V and / time series (with some U data) and the 
Meylan data provides F300W images in the first two epochs 
and F300W and F555W images in the third epoch (with limited 
time series information in each epoch). 

2.1. Astrometry 

Using the zeropoint positional offsets between the Chandra 
and //^T coordinate frames, the region within ~2f'0 of the nom- 
inal X5 position lies outside the FOV of the Gilliland data set, 
but is found on the inner part (with respect to cluster center) of 
the WF4 chip in the Meylan data. Since no Chandra source in 
the WF4 FOV with > 50 counts (not including X5) currently 
has a plausible optical counterpart, we used the PC astrometry 
to align the X-ray and optical coordinate frames (incurring a 
systematic chip-to-chip error, assumed to be Of'05, which dom- 
inates the total error budget). After this correction we found that 
only three stars are within O'.'S of the nominal X5 position, with 
separations of Of'033 (0.6a; CI), Of'23 (4.5cr; C2) and 0'/292 
(5.7cr; C3). The finding chart shows the F300W (Fig. la) and 
F555W (Fig. lb) images, from Meylan epoch 3, for the region 
around X5 and the insets in Fig. la show epochs 1 ('U(l)') 
and 2 ('U(2)'). Since faint red MS stars appear brighter in the 

STScI, which is operated by AURA, Inc. under NASA contract NAS 5-26555. 
138; pedmonds@cfa.harvard.edu; cheinke@cfa.harvard.edu, josh@cfa.harvard.edu 
I; gillil@stsci.edu 



2 



Edmonds et al. 



F555W image than in the F300W image, Figures la and lb 
show that CI has a blue color. However, C3 (just outside the 
5(7 error circle to the NE) has an even stronger blue color and 
so is also a potentially viable optical counterpart if the astro- 
metric shift between the PC and WF4 chips is much larger than 
assumed. This ambiguity is resolved by noting that CI, unlike 
C3, is clearly brighter in epoch 1 than in epoch 2 (see inset) 
confirming it as the optical counterpart (hereafter X5opt). 




Fig. 1. — Finding charts for X5 and X7. For X5 the 3rd epoch combined 
F300W ('U'; Fig. la) and F555W ('V; Fig. lb) images are shown. The X5 
companion (X5opt) is near the center of the 5a error circle, and the inset shows 
the area around X5opt from the 1st and 2nd epochs (note the clear variability 
of X5opt). For X7 the deep, oversampled combined F336W ('U'; Fig. Ic) and 
F555W ('V; Fig. Id) images from Gilliland et al. (2000) are shown. The star 
Nl is found near the center of the 20cr error circle. 

The region within a few arcsec of X7 falls on the PC images 
of both the Gilliland and Meylan data sets. Since 17 Chan- 
dra sources have likely optical counterparts on the Gilliland PC 
image, we have corrected for small linear terms in the resid- 
ual astrometric errors between Chandra and HST using least 
squares fitting. We computed the positional errors for X7 by 
adding the systematic errors to the random wavdetect er- 
rors in quadrature, resulting in la errors of 0.''0065 in RA and 
0.''0088 in Dec (see Fig. Ic and Id, where 20a error circles 
are shown). The three nearest stars to X7 in the Gilliland HST 
image are 0.''019 (2.3cr), 0.''12 (14.7cr) and 0.''23 (28.6cr) away 
(Nl, N2 and N3 respectively). Clearly, astrometrically, only Nl 
(V = 20.25; U-V= 1.72; My = 6.8) is a viable candidate for 
the optical companion of X7. Given the FOV of the PC and 
the detected number of stars on the PC chip with V < 20.25 
(6367), only 6.3 x 10"^ stars are expected within 0.''019 of Nl, 
assuming constant density over the PC FOV. 

2.2. Photometry 

The Gilliland dataset photometry (containing only X7) is de- 
scribed in GilHland et al. (2000) and Albrow et al. (2001). The 
photometry for the Meylan observations (containing X5 and 



X7) was based on combining the images at each epoch using 
drizzle routines (Hook, Pirzkal, & Fruchter 1999) in STSDAS 
and then using PSF-fitting in DAOPHOT to calculate instrumen- 
tal magnitudes. The F555W filter is a good approximation to 
Johnson V (Holtzman et al. 1995), but F300W differs sig- 
nificantly from the nearest Johnson filter (U). Therefore, we 
used ground-based photometry of 47 Tuc (Sills et al. 2000) and 
matching of main sequence (MS) turnoffs between the HST Sind 
ground-based datasets to calculate the zeropoint and then ap- 
plied corrections to F300W-y (by measuring MS ridgelines) to 
convert it toU-V. By definition this MS-ridgeline technique is 
only applied to stars (like X5opt) with colors ranging from the 
main sequence turn-off to the detected end of the MS. 




Fig. 2.— U\sU-V CMDs for the Meylan WF4 chip (containing X5opt; Fig. 
2a) and the Gilliland PC chip (containing X7; Fig. 2b). For X5opt we show 
ranges in magnitude and color given limits on the variability. Near neighbors 
to X5opt and X7 are shown (apart from C3 which is not detected in the F555W 
image. Also shown are CO WD cooling sequences from Bergeron, Wese- 
mael & Beauchamp (1995; thick lines) and He WD cooling sequences from 
Serenelli et al. (2001; thin lines). The dot-dashed sequence is the lowest mass 
model (the MSP counterpart Uopt from EGHOl is also shown). 

Since this technique is non-standard we have performed two 
consistency checks with other calibration methods. We applied 
this technique to the Meylan PC data and performed a star-by- 
star comparison between our photometry and the Gilliland et 
al. (2000) photometry. Mean differences between the two pho- 
tometric systems were < 0.05 mag in both U and V. A star- 
by- star comparison between the MS-ridgeline V calibration for 
the Meylan WF4 chip (containing X5opt) and the standard cal- 
ibration of Holtzman et al. (1995) applied to a 47 Tuc F555W 
image from the archive (program GO-6095), also gave mean er- 
rors < 0.05 mag. Combined with the 0.1 mag rms internal error 
inU-V Sit the U mag of X5opt, we estimate absolute errors for 
X5opt of ~0.2 mag in both U and V. 

The color magnitude diagrams (CMDs) for the Gilliland PC 
and Meylan WF4 images are shown in Fig. 2. Fig. 2a shows the 
mean epoch 3 CMD position of X5opt (V = 21J;U-V = 0.9; 
My = 8.2), along with reasonable ranges in magnitude and color 
given the variability (see below). With My= 8.2, X5opt has 
a similar absolute magnitude to that of the qLMXBs Cen-X4 
(My= 7.5 - 8.5; Chevalier et al. 1989) and Aql X-1 (My= 8.1; 
Chevalier et al. 1999). Also shown are CO and He WD cool- 
ing sequences (see caption). Clearly, X5opt is unlikely to be 
either a CO WD or a He WD, unlike the MSP companion Uopt 



HST detection of qLMXB companion 
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(EGHOl). Instead, it is more likely that the CMD position of 
X5opt represents the sum of a red MS star and a blue component 
from an accretion disk (see below). 
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Fig. 3. — F300W time series for X5opt for epochs 1, 2 and 3 (Figures 3a, 
b and c respectively) with time offsets between epochs as shown. For X5opt a 
simple ellipsoidal model ('ellip.'), eclipse not included, is shown (zeropoints 
adjusted to the data). Also shown is the eclipsed portion of the X-ray phase 
plot (units converted into time) from HGLOl. The \a error bars are the median 
time series rms for stars within 0.25 mag of a given magnitude level. Figure 3d 
shows the V-band time series for Nl, the nearest possible counterpart to X7. 

Fig. 2b shows the Gilliland CMD for X7 showing the posi- 
tion of the three nearest stars (Nl, N2 and N3). The only vi- 
able counterpart astrometrically, Nl, falls very close to the MS 
ridge line (also in V vs V - /) and therefore appears like a nor- 
mal MS star (unlike X5opt). This CMD position is consistent, 
within the errors, with the position in the Meylan data, and the 
F300W magnitudes from the three epochs are consistent with 
non- variability, again unlike XSopt- This suggests that Nl is 
probably not the X7 counterpart, despite the good astrometric 
match. Assuming that the real counterpart falls in the less con- 
fused part of a 5-cr error circle we set limits on its detection of 
C/ > 23, y > 23, / > 22 (using the Gilliland data). 

2.3. Time Series 

Since X5opt is near (or beyond) the limit of detectability in 
individual F300W exposures, the time series for X5opt were 
calculated by co-adding groups of 3-4 images. Figures 3a, b 
and c show the F300W time series for the 3 different epochs. 
Also shown are the eclipsed portion of the X-ray phase plot 
from HGLOl (units converted into time) and a 4.333 hour pe- 
riod sinusoid, as appropriate for X5 but with the 8.666 hour 
X-ray period divided by two to simulate a double-peaked (el- 
lipsoidal) time series. Eclipses are not included in this model. 
This model has been shifted in time and magnitude so that it 
plausibly matches the data for each epoch (the period is not 
known with sufficient accuracy to phase correct from Chan- 
dra to different HST epochs). Significant variability is seen 
within all three epochs and X5opt is clearly brighter in epoch 
1 than in the 2nd and 3rd epochs (see Fig la). This longterm 



variability is further evidence for the presence of an accretion 
disk. Note the deep eclipse observed in both epochs 2 and 3. 
The F300W eclipse appears to be significantly wider than in 
X-rays, though we do not observe the system coming out of 
eclipse. The turnover at the beginning of the second epoch 
(near Time= hr) may represent a minimum from ellipsoidal 
variability, since the timescale of variability agrees well with 
the ellipsoidal model. The 1st epoch observations may also rep- 
resent ellipsoidal variations rather than an eclipse, since in its 
brighter state the relative brightness of the disk compared to the 
secondary should be enhanced and the eclipse depth should in- 
crease. The V-band variations, not plotted here, show neither an 
eclipse nor clear evidence for ellipsoidal variations (expected to 
have a smaller amplitude than in F300W). 

No suggestion of variability is present in the Nl time series 
(Fig. 3d) and no significant signal is seen in the Nl power spec- 
trum (y or /), including the possible 5.5 hr period noted by 
HGLOl. The highest peak in the V-band corresponds to a pe- 
riod of 2.52 hours (or twice this), with a false-alarm probability 
= 0.27. The corresponding V amplitude is 0.0043 ± 0.0011 (< 
4-(7, insignificant for a blind search; similar results hold for /). 
If Nl is the X7 counterpart and is close to filling its Roche lobe 
then an inclination < 2.5° is required to reduce the amplitude 
for ellipsoidal variations from the maximum expected value of 
-^0.1, for 90° inclination, to < 0.0043. This implies that Nl is 
unlikely to be the X7 companion, and we have calculated the 
brightness limit for a faint variable star (lying near the line of 
sight of Nl) to be missed by our variability search. The X7 co- 
ordinates are so close to Nl that it will be included in any time 
series extraction. A star at V = 22.9 with intrinsic variations of 
0.1 mag superimposed on the time series of Nl would yield an 
8-(7 detection (versus the highest detected peak at ^ A-a). This 
time series limit of V ~ 23 for a companion to X7 will decrease 
for inclinations < 90°. 



3. DISCUSSION 

Using the stellar models of Bergbusch & Vandenberg (1992), 
we estimated the brightest possible secondary consistent with 
our photometry (Teff= 4100K, V = 21.7 and mass = O.53M0). 
Using the secondary radius, the X-ray luminosity of X5 and the 
binary separation (from Kepler's Third Law) we estimate that 
the maximum luminosity from heating of the secondary by the 
NS (when measured as a fraction of the secondary luminosity) 
is 2.7%. Therefore, secondary heating probably makes only a 
small contribution to the variability described above. The dom- 
inant sources of short-term variability are likely to be a com- 
bination of eclipses of the disk and hot spot by the MS star, 
ellipsoidal variations and flickering. Further observations are 
required to better define this variability. 

The likely presence of an accretion disk in X5opt, from vari- 
ability and the blue color, may appear to be inconsistent with 
the lack of X-ray evidence for accretion currently in the X5 
system, which should yield either long-term X-ray variations 
or a power law component, neither of which are seen (HGLOl). 
One possible explanation is that the X5 secondary is no longer 
filling its Roche lobe causing it to be detached from the disk. 
Such a disk would no longer be accreting matter from the sec- 
ondary, possibly causing it to enter a long-term quiescent phase 
with low density and little or no accretion onto the NS. The X5 
disk does appear to be relatively faint compared to 47 Tuc CVs, 
since the t/ - V color of X5opt (0.9) is much redder than that of 
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the 47 Tuc CVs VI, V2, Wl and W2 with ^/-V colors ranging 
from -1.25 to -0.4 (Edmonds et al. 2001, in preparation). 

To test this 'detached disk' theory we have estimated the 
degree to which X5opt fills its Roche-lobe as defined by F, 
the ratio between the stellar radius and the Roche lobe ra- 
dius. Using the Roche-lobe formula from Paczyhski (1971) 
{r/a = 0.462[(Mopt/(MNs + Mopt)]^/^ where r is the Roche-lobe 
radius, a is the binary separation, Mopt is the mass of X5opt, 
and Mns = 1.4M0), the stellar radius for a 4100 K model and 
the binary separation, we find that X5opt has F = 0.6, under- 
filling its Roche lobe. Fainter cooler secondaries will underfill 
their Roche lobes by slightly larger amounts (e.g. a star with 
My = 10.0 has F = 0.5). This behavior is consistent with the 
'detached disk' theory given above, but would be inconsistent 
with the possible detection of ellipsoidal variations of relatively 
large amplitude, requiring the secondary to have F r^l.O. The 
latter possibility would suggest that the X5 secondary is either 
bloated or slightly evolved, as appears to be the case for some 
of the CVs in NGC 6397 (Grindlay et al. 2001, in preparation) 
and as might be expected for a star undergoing mass loss. 

If the 5.5 hr X-ray period for X7 (HGLOl) is real, and if 
the X7 secondary underfills its Roche lobe by about the same 
amount as X5opt, then My~ 10.6 and V ~ 24.1, beyond our 
variability detection limits and beyond our CMD search limit 
except with the presence of a reasonably bright disk. The close 
proximity of Nl only ~0.''02 from the line of sight of X7, 
clearly makes prospects for such an optical identification dif- 



ficult. Alternatively, the period could be longer and Nl could 
be the counterpart, however the complete lack of evidence for a 
disk or variability rules against this possibility. 

The logical follow-up to these observations are spectroscopic 
studies to measure the radial velocity amplitude of absorption 
lines from X5opt. This, combined with the known inclina- 
tion and spectroscopic and photometric determinations of Mopt 
would give an estimate of the mass of the X5 NS. Using the X- 
ray spectrum constraints on the NS radius and redshift (HGLOl) 
combined with the NS mass would give the first compelling test 
of the equation of state of a NS. Also of interest would be the 
detection of emission lines in the optical spectrum. Since there 
is evidence for an accretion disk (from this work) and hot gas in 
the system (from the X-ray light curve), we expect strong disk 
or coronal emission lines to be superimposed on the absorp- 
tion line spectrum of the secondary. Study of the emission line 
profiles can test for mass outflow (visible as P Cygni profiles) 
from the system, or detect evidence for a bipolar jet (visible as 
broadened emission lines). 
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TABLE 1 

Positional, photometric and time series information for XSopt 



qLMXB 


RA^ 


Dec^ 


Xb 








period^ 




(J2000) 


(J2000) 










(hours) 


X5 


00 24 00.991(1) 


-72 04 53.202(7) 


60.7 


263.0 


22.5(2) 


21.7(2) 


8.67 


X7 


00 24 03.528(1) 


-72 04 51.938(6) 


498.5 


775.5 






5.50 



^Coordinates in MSP astrometric system; Freire et al. (2001) 

^ Pixel coordinates of X-ray sources using the STSDAS task METRIC applied to archival images 
u2ty0201t (X5) and u5jm070cr (X7) 

^ Epoch 3 mean magnitudes for X5opt 

^Periods from HGLOl; value for X7 is from a marginal detection of variability 



